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INTRODUCTION 


Ircirai  maintenance  :or  .ow-Iartfi  crcicma  LEO  < 
satellites  can  be  accompli  shed  using  forced  Kepiarian  notion 
:r  by  cang-cang  ccntroi.  Jorcea  Kepienan  notion  utilizes 
oonstant  tnrust ,  with  magnitude  equal  to  drag ,  to  prevent 
orbital  decay  and  maintain  orbital  altitude .  -Bang-bang 
control  allows  the  satellite  orbit  to  decay  a  no  -.ten  :  ires 
thrusters  to  rsooost  the  satellite  ana  keep  it  within  - 
prescribed  radial  band.  This  can  be  modeled  as  a  series  or 
orbital  transfers  offsetting  decay  due  to  atmospheric  drag. 
Currently,  LEO  satellites  are  not  required  to  maintain  a 
specified  orbital  band,  and  are  reboosted  only  to  prevent 
reentry  or  to  change  orbit. 

Historically,  co-planar  orbital  transfers  have  been 
optimally  accomplished  using  two  or  three- impulse  Hohmann 
maneuvers  ]Ref.  l:pp  78-88].  These  maneuvers  are  designed 
for  one-time  transfers  between  exo-atmospher ic  orbits,  not 
for  a  series  of  maneuvers  to  offset  drag  and  maintain  an 
orbital  band.  Bang-bang  control  utilizes  one  thruster 
firing  per  reboost  of  the  satellite,  vice  two  or  three 
firings  used  in  a  Hohmann  transfer. 

In  Ref.  2,  the  problem  of  minimum-fuel  orbital 
maintenance  for  low-altitude,  non-lifting  bodies  was 
considered,  with  the  conclusion  that  fixed-angle,  transverse 


bana-cana  control  with  the 


thrusting  i.e. 


thrust  vector 


i  constant  angle  relative  to  the  local  horizontal ,  _c  .-vs 
efficient  than  forced  Kepierian  notion.  This  tonic  ..as 
: irst  explored  bv  Ross  ana  Melton  'Ref.  ;  ]  vnc  concluded 
that  forced  Kepierian  notion  is  not  the  optimal  solution  tor 
this  problem.  Development  of  long-term,  low-ait  itucie 
satellites,  sucn  as  the  space  station,  dictate  continued 
research  into  this  topic  to  resolve  this  apparent 
contradiction. 

Little  research  has  been  done  on  the  optimisation  of 
bang-bang  control  in  relation  to  atmospheric  effects  on  non- 
lifting  (blunt}  bodies.  A  better  understanding  of  these 
effects  is  essential  for  optimization  of  orbits  and 
propulsion  equipment  for  long-term  satellites,  including  the 
space  station.  Additionally,  design  and  operation  of 
lifting  bodies  may  benefit  from  this  study. 
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II . 


GENERAL  FORMULATION  OF  THE  PROBLEM 


T 

wo  me the as 

or  optimizing  orbital 

t runs  r  e 

a  tSl 

lites  nave 

teen  studied  in  Rer._ 

and  Re i 

3  3  JL  S 

ci  numeric 

ai  simulations,  Pauis 

Rer  .  I 

"hat  banq-eang  control  is  less  erficient  than  forced 
Keplerian  notion,  and  that  transverse  thrusting  at  a  rixea 
angle  of  70  degrees  is  generally  the  optimal  feasible 
solution  far  orbital  maintenance  if  bang-oang  control  . s 
used.  This  apparently  contradicts  the  analytical  results 
Ross  and  Melton  ;  Ref .  3],  who  conclude  that  forced  i'.epieriar. 
notion  is  not  optimal  and  that  a  bang-bang  solution  exists 
that  optimizes  fuel  usage.  This  leads  Pauls  to  conclude 
that  variable  thrust-vectoring  is  essential  for  optimality 
during  finite-burns.  However,  Pauls  study  was  limited  to  a 
small  class  of  satellites  (essentially  space  platforms j  :.n 
the  upper  atmosphere.  Therefore,  ms  conclusions  are  not 
globally  valid,  and  questions  remain.  Is  transverse 
thrusting  less  efficient  than  forced  Keplerian  motion  for 
aii  parameters?  What  exactly  are  these  parameters?  Is  the 
70  degree  thrust  angle  always  optimal?  This  thesis  will  try 
to  answer  some  of  these  questions. 

The  problem  is  defined  as  maintaining  a  spacecraft 
within  a  prescribed  radial  band  R^^d. R< R^ •  Comparison  of 
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tuei  u cage  for  cana-bang  control  ana  forced  r'epier :  ar. 
.".or  ion .  over  a  wide  range  of  parameters,  1  _  1  rev:  a 
insights  on  the  optimality  question.  -.iso ,  .:eterr,:r.  : 


a  thrust  angle,  or  range  of  angles,  that  vpt  ;m.i  res  r  r.v  . 
rang  control  solution  is  aesirea. 

Ross  and  Melton  Ref.  1'  ease  their  tone i us i ors 
methods  utilizing  optimal  control  theory,  wnich  require 
initial  and  boundary  conditions  to  facilitate  a  solution 
using  numerical  integration.  This  is  a  cuiaoersor.e  ooiut 


technique  that  requires  reccmputat ion  tor  iiiar.t  vnanaes  .  r. 
initial  or  boundary  conditions.  While  this  may  te  the  most 
accurate  solution  technique,  it  is  not  easy  or  verm 
practical . 

Pauls  '.Ref.  2]  develops  a  computer  modei  to  soive  this 
problem.  Equations  of  motion  were  derived  and  solved  for 
various  parameter  combinations  to  optimize  bang-uang 
control.  Attempts  to  define  a  radial  band  and  then  control 
the  spacecraft  within  that  band  were  unsuccessful.  However, 
maintaining  specific  energy  allowed  a  band  to  be  maintained, 
albeit  larger  than  the  desired  band.  All  solutions  to  this 
model  predicted  that  forced  Keplerian  motion  is  more 
efficient  than  bang-bang  control. 

Further  study  using  a  computer  model  is  dictated.  Non- 
dimensionalization  of  the  equations  of  motion  will  allow 
effects  of  variations  in  parameters  to  be  more  easily 
analyzed.  Variations  of  parameters  will  be  analyzed  for 


» 


» 


4 


Ill 


FORMULA  DEVELOPMENT 


A.  DEVELOPMENT  OF  THE  EQUATIONS  OF  MOTION 

For  simplification,  orbital  motion  is  assumed  to 
uopianar  and  initially  circular.  Assuming  a  non- !  nt. 

riunti  body,  drag  becomes  the  only  aerodynamic  force 
perturmng  the  orbit  and  thrust  is  the  external  rorce 
applied  to  counter-balance  drag.  The  external  forces 
or.  the  spacecraft  consist  of  gravity,  aerodynamic  urag 
thrust.  Figure  l  snows  the  coordinate  system  and  net 
-  or  this  two-boav  problem. 


Transverse  axis 


Earth 


Figure  l  Graphical  Representation  of  Coordinate  System 


The  eauaticns  oi  motion  :  o; 


,v:o-:;oav 


;  l) 


(2) 


v;nere  a.  and  a.,  are  the  radial  ana  transverse  accelerations 
and  I  F.  and  r  F.r  are  the  suns  of  the  external  rad  12-  ana 
transverse  forces,  respectively,  and  n  is  spacecraft  mass. 
From  Figure  1,  the  components  of  drag  are 

(3) 


-  --sm 


/ 1 


..  -  -Dc os  ( ) 


and  the  components  of  thrust  are 


T,.=Tsin  (a) 


(4) 


(5) 


T.r=Tcos(a)  ( t 

The  angle  7  is  the  flight  path  angle  and  is  defined  cy  the 
intersection  of  the  velocity  vector  and  transverse  axis, 
while  angle  a  is  the  thrust  angle  and  is  defined  by  the 
intersection  of  the  thrust  vector  and  transverse  axis.  The 
equations  of  motion  can  now  be  written 


,  u  0.  T, 

r''  -f)  1  -r~-  —  *_ 

r~i  m  m 


(7) 


D  T 

0"r+20'  r'  =  — 

m  m 


(8) 
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where  u  is  the  Earth's  gravitational  constant.  •<  _s  the 
angular  position  or  the  spacecraft,  ana  prices  je-nczo 
differentiation  with  respect  to  tine. 

B.  NONPIMENSIONALIZATION  OF  THE  EQUATIONS  OF  MOTION 

The  equations  or  notion  are  ncnaimensiona — :oa  ..  o  that 
variations  in  parameters  are  minima  1 i zed .  The 
nonaimensionaliz ing  constants  are  chosen  to  "baiance"  the 
equations  and  thus  produce  the  same  order  of  magnitude 
responses.  This  facilitates  study  of  the  effects  or 
changing  parameters  on  the  optimization  prodiem. 

1.  Definitions 

In  order  to  nondinensionaiize  the  equations, 
variables  must  be  defined  and  nondimensionalizmg  constants 
must  be  chosen.  Nondimensionalizing  constants  for  radius 
and  time  are  designated  R  and  r,  respectively,  and  were 
chosen  somewhat  arbitrarily  based  on  the  Earth's  radius. 
Thus,  R  is  the  radius  of  the  Earth  and  r  is  the  period  of  a 
circular  orbit  whose  radius  equals  R,  and  is  defined  as 

r-2n\ir  <9> 

These  constants  define  the  nondimensionalized  variables  for 
radius  and  time 
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UO) 


rhe  nonainvens icna i is mg  factor  for  .".ass  _s  .'■!  spacecraft 
.mriai  ".ass;  ,  which  aefir.es  the  nonri mens lona^.cea  .".ass 


The  variable  >'  is  already  dimensionless. 

2.  Nondimensionalization 

Using  rhe  above  definitions,  the  following 
relaticnshiDS  are  develODed. 


dr  _  t  dr 
~d£  Rdt 

d2r  _  d  (  t  dr .  _  jn  d2r 
~dt2  ~dZ  '  ~R  dt  R  dt: 


_d£  =T£0 

dt  dt 


d2e  =  d  f  d$)  _r2  d26 

dt2  dt  dt  dt2 


Substituting  these  equations  into  the  equations  7  and  8 
substituting  equation  9  for  r  where  convenient,  and 
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rearranging  terms,  yields  the  following  nonaimensxcnai izea 
results 

d:r  -in-  _  r  Dsin-y  ,  ■ :  Tsir.a  (17 

dT-  '  ' '  dt'  r-  *  xM  * 


d2fl  =_  2  aOdr  _  r:  Dcosy  ^  r:  Tcosa  (18) 

dF  ?  dF  Hr  mM  Hr  mM 

C.  NONDIMENSIONALIZATION  OF  ATMOSPHERIC  DRAG 

1.  Definitions 

Drag  is  given  by 

D*^sre,cy-  (19) 

where  p  is  atmospheric  density  (altitude  dependent)  ,  Sre,  is 
the  reference  surface  area  of  the  spacecraft  affected  by 
atmospheric  density,  Cd  is  the  coefficient  of  drag,  and  v  is 
spacecraft  velocity. 

2.  Nondimensionalization 
Nondimensionalization  of  v2  is  accomplished  using 

the  relationship 

v2=r,2+r20,2=/ ^£)2+r2(^  )2  (20) 

1  dt  dtJ 

Substituting  equations  10,13,  and  15  gives 
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..'men  i. 3  then  suDsti'utea  into  equation  13  ana  iea 

-nto  the  torn  used  m  equations  17  and  13  to  give 

t*  D  _  7~  V  _  fjRV "  (22) 

^  2  SMt-  2EB 

where  B  is  the  initial  ballistic  coefficient  and  defined  as 


This  result  is  then  put  back  into  equations  17  ana  13  to 
yield  the  nonaimensionalized  equations 

d:r  _-^d0  4 ti"  _  pRv' siny  +  r:  Tsina  ^ 

eft2  dTT  ~W  2  mB  JnM 

d29  2  d9dr  _  pRvzcosy  +  t:  Tcosa  (25) 

dt2  r  dtdt  2mBr  Rr  mM 

Inspection  of  these  equations  reveals  that  further 
simplification  is  possible.  Assuming  an  exponential 
atmospheric  density  model  (which  will  be  used  throughout  the 
study  of  this  problem)  ,  with  atmospheric  scale  height  j 
defines 


P=P0e 


-p(r-fir.f) 


This  allows  the  ballistic  coefficient,  B,  to  be 
nondimensionalized  as 
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(27) 


Pr* 

The  thrust:  terms  can  he  similarly  r.ondir.ensicna.  .tea  by 
defining  a  nonaimensicnaiizea  thrust 

(28) 

:-ir 

Substituting  equations  27  and  28  into  equations  24  ana  25, 
ana  nondimensionalizing  the  exponent  for  iraa  result  ..n  the 
final  equations 

4;c-  _  e _ "  -.-sir.y  -aims  (29) 

r2  .a 

-gflfr- _ 

d20  _  2  dQdz  _e_ _ *  v^cosy  .  .ccsa  (30) 

St1  r  dccft  2mz3  -tr 

These  equations  will  be  used  to  generate  the  computer  model 
to  study  the  optimality  problem  of  bang-bang  control.  The 
only  parameters  that  need  to  be  varied  to  study  the 
optimality  problem  are 


a ,  p ,  T.  B 

Since  the  only  variable  parameter  in  the  nondimens lonalized 
thrust  term  is  thrust,  T,  either  one  can  be  varied  for  the 
study. 
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D.  DEVELOPMENT  OF  MASS  EQUATIONS 


~ne  or  the  optmaiity  criteria  to  ls  ..cea  m  mm  .-tuuv 
_s  one  mass  or  fuel  consumed  to  maintain  cmt.  The 
reiationsnip 

m  =  ~  (31) 

-P'" 

defines  the  rate  of  fuel  consumption  and  can  ce  mtegratea 
over  time  to  determine  the  total  fuel  consumed.  Ts  mg 
equations  11,  12,  and  29,  this  can  ce  nondimens icna i ized  ana 
simplified  to 


(32) 


fr  .  -spj 

These  equations  will  be  used  to  calculate  the  fuel  used  to 
maintain  orbit  using  bang-bang  control  for  comparison  to 
fuel  used  in  forced  Keplerian  motion. 
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IV.  DEVELOPMENT  AND  VALIDATION  OF  COMPUTER  MODEL 


A.  COMPUTER  PROGRAM  DEVELOPMENT 

A  computer  program,  written  in  FORTRAN  and  lasted  _n 
Appendix  A  with  sample  input  and  output  files,  was  level opea 
to  simulate  spacecraft  orbital  motion.  A  fourth-order 
Runge-Kutta  numerical  integration  routine  is  used  to 
integrate  the  equations  of  motion.  The  program  is  comprised 
of  a  main  program  section,  which  controls  input,  output,  and 
flow  of  information,  and  five  subroutines  that  provide 
computations  necessary  to  simulate  orbital  motion  and 
calculate  the  osculating  orbital  parameters. 

The  non-dimensionalized  equations  previously  developed 
(equations  29  and  30)  are  simplified  for  use  in  the  computer 
program  using  the  following  state  variable  definitions. 


x,  =r 

(33) 

&x\  _  dr 

(34) 

at  at 

£ 

II 

(35) 

(36) 

at  at 

These  state  variables  are  nondimensionalized  as  follows 


» 
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-  X, 


(3?) 


/ r  1%  at 


(38) 


3t 


-r  3JL  -T  Y. 

'  of  dt 


(39) 


The  state  variable  x,  is  already  nondimensionalizea .  These 
equations  are  then  substituted  into  equations  29  and  30  to 
give  the  equations  of  motion  used  in  the  computer  program. 


ax, 

H€ 


■jR<r- 


=x,xi2---L 

:<r 


*  v-sinm  Tsinct 


(40) 


2  JnE 


, -j \RCr--£)  , 

ax]  2  .—-■■■  _  e  *  vcosy  +  Tcosa  (41)  a 

-=.*x*2  2mrS  mr 

The  first  subroutine,  DRAG,  calculates  the  atmospheric 
drag  experienced  by  the  spacecraft.  A  constant  atmospheric  I 

density  model  is  used  for  initial  program  validation,  after 
which  an  exponential  density  model  is  used.  The  second 
subroutine,  EQN,  updates  the  equations  of  motion  that  define 
the  spacecraft's  orbital  motion,  equations  40  and  41.  The 
third  subroutine,  ORBPAR,  calculates  the  osculating  orbital 
parameters  for  the  spacecraft's  motion,  including  semi-major 
axis,  apogee,  perigee,  and  period.  The  fourth  subroutine, 

RK4 ,  is  a  standard  fourth-order  Runge-Kutta  numerical 
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integration  routine  that  integrates  the  state  variaoies. 

The  ^ast  subroutine.  THRUST,  governs  the  activation  ana 
deactivation  of  the  thrusters  to  maintain  the  satellite's 
:rbit . 

B.  PROGRAM  VALIDATION 

1.  Initial  validation  with  no  external  forces 

Program  validation  proceeded  in  a  logical  series  of 
steps,  with  each  step  requiring  positive  validation  before 
proceeding  to  the  next  step.  Initially,  all  external  forces 
except  grav  ty  were  neglected  in  order  to  maintain  a 
circular  orbit.  As  shown  in  Appendix  B,  Figures  1-4, 
radius,  velocity,  angular  momentum,  and  specific  energy  all 
remained  constant.  Initial  conditions  of  an  elliptic  orbit 
were  then  input  into  the  program  with  all  external  forces 
except  gravity  neglected.  Figures  5-8  of  Appendix  5  show 
that  semi-major  axis  (SMA),  eccentricity,  specific  energy 
and  angular  momentum  are  constant.  This  validates  the 
initial  program. 

2 .  Validation  of  drag 

A  constant  atmospheric  density  was  then  introduced 
into  the  program,  with  a  value  that  would  produce  a 
noticeable  effect  in  a  small  number  of  orbits  to  enhance  the 
validation  process.  The  program  was  run  for  ten  orbits  and 
the  results  compared  with  analytical  results  obtained  using 
equations  from  Ref.  4, 


•j 
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.■.men  snow  the  changes  in  semi-major  axis,  velocity,  ana 
period  per  orbit  of  a  spacecraft  under  drag.  In  ail  or 
these  equations,  ballistic  coefficient  (3)  './as  eucstitutea 
:cr  :i/c.Sr,..  The  program  was  run  using  both  raaius  • na  jemi¬ 
ma  nor  axis  for  the  term  "a"  in  the  above  equations.  Results  • 

should  be  very  similar  since  semi-major  axis  ana  radius  are 
assumea  to  be  equal  in  the  formulation  of  the  analytic 
equations  and  should  be  relatively  equal  for  the  small 
changes  expected.  Percent  errors  were  calculated  for  the 
computed  vs  analytical  results  for  changes  in  semi-major 
axis  (SMA),  velocity  (VEL) ,  and  period  (PER),  as  well  as  for 
the  percent  differences  between  computed  and  analytical 
values  of  semi-major  axis,  velocity,  and  period  after  each 
orbit.  Table  I  shows  the  results  when  radius  is  used  in  the  •; 

analytical  equations  for  "a”,  while  Table  II  shows  the 

j 

results  when  semi-major  axis  is  used  in  the  analytical 

equations  for  "a".  As  shown  in  Table  I,  the  percent  error  * 

for  the  changes  in  semi-major  axis  (%ASMA  error)  and  period 
(%APER  enor)  start  small  after  one  orbit,  grow  through  the 
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Table  I  COMPUTED  AND  ANALYTIC  ERRORS  USING  RADIUS 


ORBIT 


i  iap 

;  3  PER 

.AY 
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!  ERROR 

:  ERROR 

,  ERROR 
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6  10.017 
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.3017 


. 00005 
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00011  |  11. 666  i  . 00020 
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0  00  5  6 


00069  i 


00031  1 


0009- 


00 1'0~  ' 
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Table  II  COMPUTED  AND  ANALYTIC  ERRORS  USING  SEMI-MAJOR  AXIS 


ORBIT 

/ 

%ASMA 

ERROR 

%  SMA 
ERROR 

%AP 

ERROR 

%  PER 
ERROR 

%AV 

ERROR 

%  YEL  | 

ERROR  | 

1 

9.5531 

. 00002 

7 . 2375 

.00015  ! 

2 

9.5531 

. 00003 

9.5532 

. 00004 

5.1534 

.00029  !| 

3 

9.5531 

. 00004 

9.5532 

.  00006 

4 

9.5532 

.  00005 

9.5533 

.  00007 

2 .2979 

.  00056 

5 

9.5532 

.  00006 

9.5534 

. 00009 

6.6409 

.  00069 

6 

9.5532 

.  00007 

9.5535 

. 00011 

10.484 

.  00081 

7 

9.5533 

. 00009 

9.5535 

. 00013 

13 . 155 

.  00094 

8 

9.5533 

.  00010 

9.5536 

.  00015 

14 .265 

.  00107 

9 

9.5533 

.  00011 

9.5537 

.00017 

13 .859 

.  00120 

10 

9.5534 

.00012 

9.5537 

.00019 

12 . 377 

.  00134 
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■r;jc;veiv  constant  .IT  c-rror  i  our.  a  oner.  us;na  -er.i--.g-'.  or 
axis  m  Table  TI .  However,  _n  both  ;ases ,  >s  cncwn  .r. 
rabies  I  ana  II,  the  percent  error  between  ~r.e  oaxcu iatea 
ana  analytic  semi-aaicr  axis  .  1  SMA  error;  ana  period  1 E? 
error;  after  eacn  orbit  is  less  than  0.00031.  The  absolute 
percent  error  in  the  change  m  velocity  IAV  error;  ana  the 
iifference  between  the  calculated  ana  analytic  eiccity 
"EL  error:  after  eacn  orbit  are  the  care  .-.-nether  icing  cen¬ 
tal  or  axis  or  radius  for  "a",  as  shown  in  Tables  :  ana  71. 
Again,  vnile  the  percent  error  in  the  change  in  velocity 
grows  to  greater  than  ten  percent,  the  percent  error  between 
calculated  and  analytic  velocity  after  each  orbit  very  snail 
less  than  0.002%).  The  reason  for  the  ten  percent  errors 
between  parameter  changes  appears  to  be  the  accuracy 
involved  when  the  changes  are  so  small  { i . e.  seri-maior  axis 
is  changing  by  approximately  200  meters  per  orbit  out  or  an 
initial  orbit  of  6638200  meters) .  Thus  an  error  of  10%  in 
change  of  semi-major  axis  is  only  20  meters,  which  is 
insignificant  in  terms  of  the  initial  radius.  The  analytic 
and  computed  results  are  thus  in  agreement  and  the  drag 
portion  of  the  program  is  validated. 
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3.  Validation  of  thrust 

The  -hr use  portion  or  the  program  was  vai.aat 
throughout  the  study  ov  having  drag  present  an a  toe 
spacecraft  maintain  an  orbitai  band  by  reboosting  usi 
thrusters . 

C.  ORBIT  CONTROL  STRATEGY 

It  is  desired  to  design  a  control  strategy  that  . 
simple  and  will  maintain  the  spacecraft  within  a  pre¬ 
specified  radial  band.  Figure  2  shows  a  typical  band 
defining  limits  of  radius  and  specific  energy. 


Figure  2  Orbital  band  with  radius  and  energy 


Ail  control  strategies  were  run  using  numerous  cor.Dir.at.ions 
or  parameter  values  f  i.e  thrust  vector  angle.  . ,  '..'as  "anea 
crom  0-90°  in  5'  increments,  radial  band  was  vanea  using 
15  ana  100  km  sands .  thrust  was  varied  from  2  5  to  ;00  , 

etc.. .  For  consistency,  all  cases  presented  here  were  run 
vith  the  following  baseline  values 

R0=6 6  3  8. 2km 
M-20000/rgr 
Thrusc=3Q0N 

3a. ~ lisa ccosffi ci c -3- 1 3 C kz  ~ 

pn=1i,**a^“-.wiCCr  .7? " 

j3  =2  . 12 e-zkm  • 
rsp-3  0Csec 

a  -7  0° 

Radialband=2km 

which  produce  representative  results  corresponding  to 

3=25000 
T=0 . 060 

1.  Control  using  radius 

Initial  attempts  were  made  to  control  the  radial 
band  using  radius  as  a  parameter,  i.e.  turn  thrusters  on 
when  the  radius  drops  to  the  low  end  of  the  band  (RMIN)  and 
turn  them  off  when  the  satellite  is  reboosted  to  the  high 
end  of  the  band  (RMAX) .  As  shown  in  Figures  1  of  Appendix 
C ,  this  control  strategy  was  ineffective  because  the 
specific  energy  of  the  orbit  continually  increased  and  never 
leveled  off  (see  Appendix  C,  Figure  2) .  The  eccentricity  of 
this  orbit  grew  to  approximately  0.1. 
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2.  Control  using  total  energy 


Control  was  next  attemptea  using  total  energy.  The 
total  energy  of  the  initial  orbit  was  calculated,  as  '..ere 
the  total  energies  of  the  spacecraft  in  circular  omits  at 
the  minimum  and  maximum  of  the  radial  band  (see  Figure  Z) 
using 


The  thrusters  were  turned  on  when  the  spacecraft  total 
energy  dropped  below  the  minimum  orbit  total  eneray  ana 
turned  off  when  the  spacecraft  total  energy  rose  above 
maximum  orbit  total  energy.  This  method  failed  to  control 
radius,  as  the  spacecraft  orbit  always  eventuall”  decayed, 
but  successfully  controlled  energy,  as  shown  in  Figures  3 
and  4  of  Appendix  C.  While  this  method  succeeded  in 
maintaining  a  pre-determined  radial  band,  the  orbital  radius 
was  not  maintained  and  decayed. 

3.  Control  using  radius  and  total  energy 

A  combination  control  strategy  utilizing  both  radius 
and  total  energy  was  attempted.  Thrusters  were  turned  on 
when  the  orbital  radius  dropped  to  the  minimum  radius  and 
turned  off  when  total  energy  increased  to  the  initial  total 
energy.  Figures  5  and  6  of  Appendix  C  show  the  unsuccessful 
results.  Total  energy  would  seem  to  be  the  parameter  of 
choice  for  maintaining  the  orbit  vice  specific  energy. 


»1 
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However,  if  the  total  energy  is  maintained  at  a  constant 
value  the  orbital  radius  must  decrease  and/or  velocity  must 
decrease  since  the  spacecraft  mass  decreases  due  to 
thrusting.  This  can  be  shown  using  equation  13. 

4.  Control  using  specific  energy 

The  next  attempt  at  control  involved  spacecrart 
specific  energy.  This  control  was  similar  to  the  attempt 
using  total  energy,  only  specific  energy  was  used.  Specific 
energy  is  a  function  of  orbital  radius  (for  circular 
orbits) ,  so  "conserving"  specific  energy  should  result  in  a 
"conservation"  of  radius  (i.e.  a  constant  average  radius) . 

As  shown  in  Figure  7  of  Appendix  C,  the  orbital  radius  was 
maintained,  however  the  orbital  band  was  not.  Conserving 
spacecraft  specific  energy,  as  shown  in  Figure  8  of  Appendix 
C,  failed  to  maintain  the  initial  radius  of  the  spacecraft. 

5.  Control  using  radius  and  specific  energy 

Control  using  radius  and  spacecraft  specific  energy 
is  successful  in  maintaining  radius  and  a  radial  band,  as 
shown  in  Figures  9  and  10  of  Appendix  C.  Thrusters  are 
turned  on  based  on  the  spacecraft  orbit  decaying  below  a 
specified  radial  minimum  and  spacecraft  specific  energy 
being  below  the  initial  specific  energy.  Thrusters  remain 
on  until  specific  energy  increases  above  the  initial  value. 
This  control  strategy  results  in  an  orbit  whose  radius 
oscillates  about  the  initial  radius,  but  not  necessarily  at 


the  pre-specif led  value.  The  variation  of  the  actual  raaxai 
oand  from  the  desired  band  varies  and  is  exanunea  in  greater 
detail  in  the  next  chapter. 


» 
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V.  ANALYSIS  AND  RESULTS 


Analysis  of  the  validated  control  strategy  is  performed 
in  five  steps.  A  set  of  parameter  values  is  chosen  as  the 
baseline  case  about  which  they  will  be  varied.  These  values 
are  selected  as  representative  of  a  space-station/platform 
type  satellite.  The  baseline  values  are 


R- 6 6  3  3. 2  km 
M=20QQ0kg 
Isp= 3  00  sec 
r=300w 

p0=9  .  Ze-10kg/ mz 
B=l50kg/m2 
5=25000 
p=2.l2e-5ion‘1 


The  simulation  program  is  run  for  100  orbits,  with  the 
equations  of  motion  updated  5000  times  per  orbit,  and  output 
sampled  100  times  per  orbit. 


A.  VARIATION  OF  THRUST  ANGLE,  a 

First,  numerous  cases  are  studied  to  determine  the 
validity  of  an  optimal  thrust  angle,  a,  of  70°.  The  program 
is  run  with  combinations  of  variations  of  parameters  about 
the  baseline  values,  while  varying  a.  from  0-90°  in  5 ° 
increments.  In  all  cases,  thrust  angles  of  65  and  70°  are 
the  only  values  that  maintain  a  constant  radial  band. 

Figures  1-4  of  Appendix  D  show  results  over  the  range  of  60- 
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7 5°  for  the  baseline  parameters,  and  are  representat :v 
all  cases.  Values  of  a  below  60”  and  above  "5”  produce 
results  that  degrade  significantly  from  those  sncv/n.  -sr.  -  ~ 
major  axis  (SMA)  and  eccentricity  are  also  compared 
throughout  the  analysis  for  determining  optimal  values  or  j.  . 
Figures  5  and  6  of  Appendix  D  show  results  consistent  v/itn  a 
maintained  radial  band  (a=70°)  ,  while  Figures  "  and  3  show 
results  consistent  with  a  radial  band  not  being  maintained 
(a:=50o)  .  As  can  be  seen,  eccentricity  must  ^evel  off  to  a 
constant  value  to  maintain  a  radial  band,  v hi.e  SMA  aicne  ~s 
not  an  indicator  of  orbital  maintenance. 

The  mass  of  fuel  required  to  maintain  orbit  generally 
increases  as  a  increases  from  0  to  90°,  with  some  minor 
local  fluctuations.  Table  III  shows  the  mass  of  fuel  burned 
in  kilograms  (kg)  to  maintain  2,  25,  and  100  km  pre¬ 
specified  bands,  with  a  varying  from  60-75°  and  baseline 

Table  III  MASS  OF  FUEL  BURNED  (kg)  _ ___ 


Band  (km) 

a=60° 

a=65° 

a=7  0° 

a=7  5” 

2 

2049 

2013 

2174 

2683 

25 

2643 

2828 

3243 

3595 

100 

3020 

2853 

2826 

3787 

values  for  other  parameters,  and  is  typical  of  the  trend 
seen.  Thus,  from  a  fuel  usage  standpoint  a  should  be 


minimized,  but  from  a  radial  band  standpoint  a  must  be  in 
the  range  of  65-7  0°. 

This  shows  that  a  constant  thrust  angle,  a,  or 
approximately  65-70°  is  necessary  to  maintain  a  radial  band. 
The  following  analyses  are  conducted  with  a  range  of  60-75’ 
for  a  to  minimize  the  data  collection  while  ensuring  the 
continuing  validity  of  this  result. 

B.  VARIATION  OF  PRE-SPECIFIED  RADIAL  BANDWIDTH 

Once  a  is  chosen,  it  is  desired  to  know  how  the  pre¬ 
specified  radial  bandwidth  affects  the  results.  While 
maintaining  a  within  the  60-75°  range  and  using  the  baseline 
values  for  other  parameters,  the  program  is  run  using 
specified  radial  bandwidths  of  1,2,5,10,25,100,  and  200  km. 
Although  the  large  bandwidths  are  not  practical,  they  are 
used  to  determine  the  existance  of  any  limiting  conditions. 
These  results,  shown  in  Table  IV  (first  entry  is  fuel  burned 
in  kg,  second  entry  is  actual  radial  band  maintained  in  km) , 
show  that  fuel  required  to  maintain  the  orbit  generally 
increases  as  the  controlled  radial  bandwidth  increases,  with 
the  exception  being  the  100  km  band.  Although  the  error 
between  the  radial  band  maintained  and  the  pre-specif ied 
band  is  minimized  as  a  increases,  the  maintained  band  cannot 
be  correlated  to  the  pre-specif ied  band  (changing  any 
parameter  changes  the  maintained  band) .  For  this  reason, 
the  radial  bands  discussed  throughout  this  analysis  are  pre- 
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Table  IV  MASS  OF  FUEL  BURNED  fkg)  AND  ACTUAL  RADIAL  SAND 
MAINTAINED  :km) 


o 

' V 

II 

in 

II 

1537/0.7 

1684/0.7 

2049/2 

2013/1.8 

2516/12 

2628/8 

2392/25 

2554/16 

2643/55 

2828/43 

2812/90 

3049/75 

3020/130 

2859/125 

2993/130 

3549/140 

2  =  70 


2025/1 . 0 


3191/13 


3243/35 


3136/65 


2527/ 1 . 1 


2683/2 


065/5 


3602/11 


3595/30 


7  8  7/75 


i  '  t 


specified  bands.  Figures  1-3  of  Appendix  E  show  the  mass  of 

fuel  burned  using  bang-bang  control  and  forced  Keplerian  • 

control  for  2,  25,  and  100  km  bands  with  a=70°.  As  can  be 

seen,  the  fuel  used  increases  in  a  generally  linear  manner 

as  the  number  of  orbits  increases.  For  the  100  km  band,  the  V 

fuel  used  is  linear  when  the  points  at  the  end  of  each  burn 

are  connected,  as  indicated  by  the  dotted  line  in  Figure  3. 

Thus,  these  results  can  be  extrapolated  for  any  number  of  • 

orbits  to  give  fuel  usage  vs  Keplerian  usage.  However,  the 

thrusters  do  not  fire  at  the  same  point  in  each  orbit  (the 

thrusters  fire  approximately  every  three-fourths  of  an  orbit  * 

for  the  cases  of  a  small  radial  band  and,  thus  firing 

precesses  each  orbit) .  These  results  are  consistent 

throughout  the  analysis.  * 
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A  constant:  radial  band  is  maintained  for  a  cr 
approximately  65-7 cr’  but  r.ot  for  60  ana  7  5'*,  as  snown  _n 
Figures  4-7  of  Appendix  E ,  which  again  verifies  the  required 
thrust  angle.  These  figures  are  for  a  pre-specif isa  raaiai 
bandwidth  cf  25  km,  the  results  are  similar  for  candwidths 
of  1-50  km.  Figure  8  of  Appendix  E  shows  the  results  for  a 
pre-specif ied  bandwidth  of  100  km  and  an  a  of  70''  ;this  is 
representative  of  results  for  a=60-75°  ana  bands  of  100  km 
and  greater) .  As  can  be  seen,  the  band  is  wide  enough  that 
thrusters  fire  infrequently  to  reboost  the  satellite.  The 
actual  radial  bandwidth  maintained  was  then  compared  to  the 
specified  band.  Based  on  these  results,  all  subsequent 
analysis  is  performed  for  three  cases  of  radial  bandwidth; 

2,  25,  and  100  km.  This  will  give  results  for  a  small, 
medium  and  large  bandwidths,  while  examining  the  unusual 
results  at  100  km. 

C.  VARIATION  OF  ATMOSPHERIC  SCALE  HEIGHT,  0 

Of  the  three  variable  parameters  in  the  non- 
dimensionalized  equations,  0  should  have  the  most 
predictable  effect  on  the  results  and  is  analyzed  first. 
Using  baseline  values  for  other  parameters,  0  is  varied  from 
1.0e-10  to  l.0e-2,  and  results  are  examined  for  logical 
trends.  Logic  would  dictate  that  increasing  the  value  of  0 
will  increase  the  rate  of  orbital  decay  within  the  band, 
requiring  more  thruster  firings  to  reboost  the  satellite, 


and  thus  require  more  fuel.  This  is  indeed  the  case  as  me 
results  in  Appendix  F  snow.  Figure  1  shows  the  fuei  used  as 
a  function  of  3  for  a  pre-specif led  2  km  bana,  .mile  Figures 
2  and  3  show  the  same  results  for  25  and  100  km  bands, 
respectively.  All  three  plots  are  similar  in  appearance 
indicating  changes  in  3  affect  orbital  decay  uniformly  over 
the  range  of  bands.  Below  a  certain  value,  ie-5  for  2  km 
band,  changes  in  3  have  no  effect  on  fuel  usage.  Above  a 
certain  value,  le-3  for  2  km  band,  the  orbit  decays  too 
quickly  for  thrusters  to  control.  In  between  these  values, 
increases  in  3  increasingly  affect  the  fuel  used.  Since 
changes  in  3  do  not  affect  the  forced  Keplerian  solution, 
the  mass  of  fuel  required  for  forced  Keplerian  maintenance 
for  all  cases  is  approximately  649  kg.  Thus,  smaller  values 
of  3  allow  the  bang-bang  control  to  be  closer  in  efficiency 
to  forced  Keplerian  motion.  Analysis  of  the  plots  for 
radial  band  again  show  that  a  radial  band  is  not  maintained 
for  60  or  7  5°  thrust  angle,  but  is  maintained  and  optimized 
in  the  65-70°  range.  From  these  results,  it  is  shown  that 
changes  in  3  affect  the  efficiency  of  bang-bang  control  as 
expected  and  are  predictable. 

D.  VARIATION  OF  THRUST ,  T 

The  next  variable  to  be  analyzed  is  thrust  (T) .  As 
explained  in  Chapter  3,  either  thrust  or  nondimensionalized 
thrust  can  be  varied  and  the  results  will  be  the  same. 
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Thrust  is  cnosen  cecause  variations  _n  thrust  xre.  tore 
easily  related  to  ana  understood.  "or  this  analysis,  thrust 
is  varied  ever  the  range  10-1000  ,  mile  ail  other 

parameters  are  fixea  at  the  oaseiine  value.  Figures  l-i  or 
Appendix  G  show  the  fuel  required,  for  a-7 0°.  to  maintain  ~ , 
25,  and  100  km  pre-specif ied  bands,  respectively. 

Trends  are  similar  on  ail  three  graphs.  Fuel  usage 
starts  off  low  for  small  values  of  thrust,  with  thrust  being 
small  enough  to  produce  quasi-forced  Keplerian  motion 
(thrusters  fire  for  majority  of  time  to  reboost  satellite;. 
Forced  Keplerian  thrust  is  determined  to  be  approximately 
2.8  N  for  the  reference  conditions.  To  obtain  this  value  cf 
transverse  thrust  using  a  thrust  angle  of  70"  would  require 
a  thrust  of  approximately  11.1  N.  Thus,  small  values  of 
thrust  may  produce  a  quasi-forced  Keplerian  motion  with 
respect  to  the  transverse  thrust  component.  This  is 
followed  by  a  period  of  increased  fuel  usage  as  thrust 
increases.  In  the  2  km  band  case,  increasing  thrust  further 
shows  a  drop  in  fuel  usage  followed  by  a  relatively  constant 
fuel  usage  as  thrust  increases  above  300  N.  The  25  km  band 
case  shows  more  fluctuation  in  fuel  usage  as  thrust 
increases,  followed  again  by  a  relatively  constant  fuel 
usage  as  thrust  increases  above  300  N.  For  the  100  km  band 
case,  increasing  thruster  size  generally  increases  fuel 
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usage  except  for  a  fairly  constant  region  in  the  range  or 

50-250  N.  * 

Thus,  other  than  a  possible  minimizing  ci  fuei  -sage  cue 
to  this  apparent  quasi-forced  Keplerian  motion,  changes  in 
thrust  appear  to  have  little  effect  on  the  optimality 
problem.  The  smaller  fuei  usage  totals  at  low  thrust  25-50 
N)  may  be  indicative  of  approaching  a  quasi-forced  Keplerian 
motion  (except  thruster  angle  is  70°  vice  0°) .  For  all  three 
cases,  as  thrust  increases,  fuel  usage  appears  to  approacn  a 

constant  value.  An  optimal  thrust  angle  of  55-70"  is  again 

.  .  • 

verified  m  this  portion  of  the  analysis. 

E.  VARIATIONS  IN  NONDIMENSIONALIZED  BALLISTIC  COEFFICIENT 

The  last  parameter  to  be  analyzed  is  nondimensionai  * 

ballistic  coefficient,  B-.  3y  nondimensionalizing , 
variations  in  B  will  encompass  variations  in  the  individual 
parameters  of  mass,  atmospheric  density,  coefficient  of 
drag,  and  reference  surface  area,  as  well  as  combinations  of 
variations  involving  these  parameters. 

Analysis  is  conducted  by  varying  B  over  the  range  4e3-  » 

2e5  and  comparing  mass  of  fuel  burned  using  bang-bang 
control  (MFB)  to  mass  of  fuel  burned  using  forced  Keplerian 
motion  (MFBK) .  This  analysis  is  conducted  for  the  three  ft 

bands  of  2,  25,  and  100  km  using  baseline  values  for  other 
parameters  while  varying  a  from  60-75°.  Figures  1  and  2  of 
Appendix  H  show  the  results  for  the  2  km  band  case,  results  * 
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maximizing  B  minimizes  fuel  usage. 

However,  this  does  nor  give  a  feel  rcr  how  erncient 
bang-oang  control  is  compared  to  rorced  Kepienan  motion. 

To  analyze  this,  the  ratio  of  MFB  to  MFBK  is  plotted  versus 
3  for  2,  25,  and  100  km  bands  and  a  of  65  ana  70".  These 
results  are  shown  in  Figures  2-3  of  Appendix  H.  .-.gain,  3  is 
varied  over  the  range  4e3-2e5.  Values  of  B  below  and  above 
this  range  (and  the  range  plotted  for  100  km  band)  cause  the 
orbit  to  decay  too  rapidly  or  too  slowly,  respectively. 
Logically,  this  is  analogous  to  increasing  the  atmospheric 
density  to  a  value  that  causes  the  orbit  to  decay  faster 
than  thrusters  can  maintain,  or  decreasing  the  atmospheric 
density  to  a  value  that  causes  the  orbit  to  decay  so  slowly 
that  the  thrusters  never  fire  in  100  orbits  (thus  giving  no 
data  since  at  least  one  firing  of  the  thrusters  m  bang-bang 
control  is  required  for  comparison  to  forced  Kepierian 
motion) . 

As  seen  in  Figures  3-8,  the  general  results  for  a  of  65 
and  70°  are  similar  for  all  three  bands.  Also,  the  results 
for  all  three  bands  show  a  roughly  constant  value  for  the 
ratio  of  MFB  to  MFBK,  with  the  2  km  band  being  most 
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efficient  compared  to  forced  Kepierian  .tot ion.  roilowea  ;.v 
one  100  Km  cana  ana  then  tr.e  15  Km  cana.  A  vaiue  or  1  ;  or 
the  ratio  would  indicate  that  the  cang-cang  control  .as 
equally  efficient  to  Kepierian  motion.  Values  .ass  than  1 
mould  indicate  that  the  cang-bang  control  is  superior  to 
Kepierian  motion  and  mould  numerically  validate  the 
theoretical  results  of  Ross  and  Melton  [Ref.  1'.  The 
results  shown  are  consistently  above  a  ratio  of  2, 
indicating  that  forced  Kepierian  motion  is  sianif leant ly 
more  efficient  than  bang-bang  control  over  ail  values  or  2. 

Throughout  this  analysis,  the  optimality  or  a  m  the 
range  of  65-70°  was  verified.  This  result  has  remained 
valid  throughout  all  of  the  analyses  conducted,  indicating 
that  this  is  a  globally  valid  result. 

» 


» 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  purpose  of  “his  thesis  was  to  det.err.ine  if  there  ire 
sny  cases  where  fixed  angle,  Pang-bang  control,  orcitai 
maintenance  approaches  or  exceeds  the  efficiency  of  forced 
Kepierian  notion  as  proposed  by  Ross  and  Melton  ;Ref.  i". 

By  nondimensionalizing  the  equations,  a  thorough  study  of 
the  effects  of  all  parameter  combinations  was  possible. 
Analysis  of  the  results  indicate  that,  in  ail  cases,  forced 
Kepierian  motion  is  superior  to  bang-bang  control.  There 
are  no  cases  where  bang-bang  control  even  remotely 
approaches  the  efficiency  of  forced  Kepierian  motion.  Thus, 
these  results  appear  to  be  globally  valid. 

The  optimality  question  for  thrust  angle  also  appears  to 
be  resolved  globally.  In  all  cases,  the  thrust  angle  is 
required  to  be  between  65  and  70°  in  order  to  maintain  a 
constant  radial  band.  Generally,  minimizing  this  angle 
improves  fuel  consumption,  so  the  angle  should  be  minimized 
as  far  as  practicable  while  ensuring  a  constant  radial  band. 

Maintenance  of  a  pre-specif ied  radial  bandwidth  is  not 
always  possible  or  predictable  for  this  control  strategy. 

The  error  between  specified  and  maintained  band  can  be 
minimized  through  proper  bandwidth  selection.  Selecting  a 
relatively  small,  2  km,  or  large,  100  km,  band  minimizes 
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the  error.  Also,  maximizing  the  thrust  anaie,  a,  generally 
minimizes  the  error,  and  thus,  a  should  be  see  as  man  as 
possible,  while  maintaining  a  constant  radial  band.  Since 
this  contradicts  the  requirement  to  minimize  a  tor  fuel 
savings,  trade-offs  must  be  considered.  Minimizing  the 
specified  bandwidth  moves  the  efficiency  of  bana-cana 
control  closer  to  that  of  forced  Keplerian  motion. 

The  effects  of  changes  in  -mospheric  scale  height,  8 , 
are  logical  and  predictable  cn  the  effects  of  the  mass  of 
-_ai  burned  in  oang-bang  control.  Since,  by  definition, 
changes  in  ,5  do  not  affect  forced  Keplerian  results,  the 
changes  in  mass  burned  in  bang-bang  control  directly 
correlate  to  changes  in  efficiency.  Thus,  the  lower  the 
value  of  £,  the  more  efficient  bang-bang  control  is  in 
relation  to  Keplerian  motion. 

The  effects  of  varying  thrus*-  are  not  so  predictable ; 
however,  general  trends  can  be  imm_ed  from  the  results. 

Low  values  of  thrust  appear  to  approach  quasi-forced 
Keplerian  motion  and  are  more  efficient.  Higher  values  of 
thrust  increase  the  mass  of  fuel  burned,  however,  a  maximum 
limit  appears  to  be  approached.  As  in  the  case  of  8 , 
changes  in  thrust  do  not  affect  the  mass  of  fuel  burned  in 
Keplerian  control,  thus  the  mass  of  fuel  burned  in  bang-bang 
control  is  indicative  of  the  efficiency. 

This  is  not  the  case  for  changes  in  B,  where  changes 
affect  the  final  Keplerian  fuel  usage.  Thus,  in  order  to 
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examine  efficiency,  the  ratio  of  aass  of  fuei  our nea  -sir.g 
uang-oang  control  to  -ass  of  fuei  burned  v/ith  .'eoierian 
notion  must  re  used.  it  v.as  initially  expected  mat  marges 
-n  5  would  have  noticeable  effects  on  efficiency  ana  'could 
show  a  trend  that  wouid  indicate  values  where  cang-bang 
tontrol  efficiency  approached  Keplerian  efficiency.  This  is 
not  the  case,  as  changes  in  B  seem  to  have  little  effect  on 
the  ratio  of  efficiencies. 

Thus,  the  superiority  of  forced  Keplerian  motion  to 
fixed  angle  bang-bang  control  is  globally  confirmed. 

However,  further  research  should  be  conducted  into  the 
solutions  proposed  by  Ross  and  Melton  [Ref. 3’.  Control 
using  variable  thrust  angle  bang-bang  control  so-called 
Primer-vectoring)  should  be  examined  for  comparison  to  these 
results.  Also,  fixed  angle  bang-bang  control,  as  described 
here,  should  be  compared  with  other  innovative  methods  of 
control,  as  well  as  with  the  "Lambert  control"  used  with  the 
shuttle,  for  possible  fuel  savings  for  space-station  type 
platforms . 

The  results  of  this  thesis  also  suggest  areas  for 
further  study.  A  method  for  maintaining  a  pre-specif ied 
band  is  highly  desirable,  but  was  not  achieved  in  this 
study.  Although  the  global  validity  of  the  thrust  angle 
being  required  to  be  between  65  and  70°  was  shown,  the 
reason  for  this  has  not  been  studied.  The  problem  of 


micro-gravity  constraints  and  its  effect  on  thruster  nrmgs 
and  control  strategies  has  also  not  been  studied,  out  v/ouid 
be  a  critical  factor  in  any  control  design.  The  ract  that 
forced  Keplerian  motion  has  been  shown  to  be  on  the  order  or 
three  times  more  efficient  than  bang-bang  control  indicates 
that  there  is  much  room  for  research  into  further 
optimization  techniques  for  bang-bang  control.  Finally,  the 
question  of  practicality  has  not  been  addressed  m  relation 
to  both  forced  Keplerian  motion  (generally  not  very 
practicable  but  very  efficient)  and  bang-bang  control  i very 
practicable  but  not  very  efficient) . 
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APPENDIX  A 


A.  PROGRAM  LISTING 


OBJECTIVE:  COMPARISON’  OF  FORCED  KEPLARIAN  TO  BANG-BANG 
CONTROL  OF  ORBITAL  ALTITUDE 


c 

C 

C 

C 


c 

C 

C 

C 

c 

c 

c 

c 

c 

c 

r* 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


VARIABLE  DEFINITIONS 

XBAR( l)=NONDIMENSIONALIZED  ORBITAL  RADIUS 

XBAR( 2 ) =NONDIMENSIONALI2ED  ORBITAL  RADIAL  VELOCITY 

XBAR ( 3 ) =THETA 

X8AR(4)=NONDIMENSIONALI2ED  ORBITAL  ANGULAR  VELOCITY 
XBDOT ( 1 ) =DERIVATIVE  OF  XBAR(l) 

XBDOT ( 2 ) =DERIVATIVE  OF  XBAR(2) 

XBDOT ( 3 ) =DERIVATIVE  OF  XBAR(3) 

XBDOT ( 4 ) ^DERIVATIVE  OF  XBAR(4) 

RO=INITIAL  ORBITAL  RADIUS 
R=ORBITAL  RADIUS 

RREF=REFERENCE  ALTITUDE  FOR  DENSITY 
D=DRAG  (N) 

DK=DRAG  KEPLERIAN 
E=SPECIFIC  ENERGY 

EO=SPECIFIC  ENERGY  AT  INITIAL  RADIUS  RO 
M=S/C  MASS  (KG) 

MO=INITIAL  S/C  MASS  (KG) 

ME=ARBITRARY  MASS  FOR  NONDIMENSIONALIZING  (20000KG) 

MBAR=NONDIMENSIONALIZED  MASS 

MF=MASS  OF  FUEL  BURNED  IN  TIME  INCREMENT 

MFK=MASS  OF  FUEL  BURNED  PER  TIME  INCREMENT  KEPLERIAN 

MFT=MASS  OF  FUEL  BURNED  TOTAL 

MFTK=MASS  OF  FUEL  BURNED  TOTAL  WITH  KEPLERIAN  MOTION 
GAMMAR=FLIGHT  PATH  ANGLE  (RAD) 

GAMMAD=FLIGHT  PATH  ANGLE  (DEG) 

TH=THRUST  (N) 

THK=KEPLERIAN  THRUST 
THMAX=BLOWDOWN  (MAX)  THRUST  (N) 

THBAR=NONDIMENSIONALIZED  THRUST  = (TH*TAU**2 ) / (MR) 
THBM=NONDIMENSIONALIZED  BLOWDOWN  (MAX)  THRUST 
ALPHAR=THRUST  ANGLE  (RAD) 

B=BALLISTIC  COEFFICIENT  (M/(CD*S)) 
BBAR=NONDIMENSIONALIZED  B  =B/(RHOO*RE) 

RHO=CALCULATED  ATMOSPHERIC  DENSITY 
RHOO=REFERENCE  ATMOSPHERIC  DENSITY 
SPI=SPECIFIC  IMPULSE 
V*VELOCITY 

VBAR=NOND IMENS IONALI ZED  VELOCITY 
G=GRAVITATIONAL  ACCELERATION 
SMA=SEMI-MAJOR  AXIS 
ECC=ECCENTRICITY 
T=TIME 

TINC= INCREMENT  OF  TIME ( STEP  SIZE) 
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TF=END  TIME 

?TI=PRINT  TIME  INTERVAL 


C  CONSTANTS 


C 

c 
c 
c 
c 

Q 

c 

c  VARIABLE  DECLARATION 

IMPLICIT  REAL* 8 ( A-H, M-Z ) 

DIMENSION  XBAR(4) , XBDOT ( 4 ) 

C 

c 

c  CONSTANT  DEFINITIONS 

?I=DATAN ( 1 . 0D+00 ) *4 . 0D+00 
MU=3 . 98601208133D+14 
G=9 . 806D+0 
RE=6 . 3782D+6 

TAU=2 . *PI*RE**1. 5 /MU* *0.5 
J=4 

ME=20000.D+0 

C 

C  INITIALIZE  FUEL  USAGE  TOTALS 

C 

MFT=0 . QD+0 
MFTK=0 . 0D+0 
C 

c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
c  MAIN  PROGRAM 

c 

OPEN { 10 , FILE= ' inp ' , STATUS= ' OLD ' ) 

OPEN ( 20 , FILE= ' outp' , STATUS='NEW' ) 

OPEN ( 30 , FILE= ' orbp ' , STATUS =' NEW ' ) 

OPEN ( 40 , FILE* ' ratp ' , STATUS =' NEW' ) 

C 

C  READ  INPUTS  AND  INITIALIZE  PARAMETERS 

C 

READ ( 10 , 1 ) R0 , M0 , THMAX , TB , TFB , TINCB , PTIB , BBAR , BETA , SPI 
1  FORMAT(10(/, 21X,D13.7) ) 

C 

C 

M=M0 

MBAR-M/ME 

TAUR0=2.*PI*R0** 1.5/MU* *0.5 

TINC=TINCB*TAU 

THBM=THMAX  *TAU  *  TAU / (ME*RE ) 

C 

C  READ  PARAMETERS  THRUST  ANGLE,  DESIRED  BAND, THRUST  FACTOR 

C 

c  PRINT*, 'ENTER  ALPHA' 


RE=RADIUS  "  EARTH 

MU=EARTHS  GRAVITAIONAL  CONSTANT 

TAU=ORBITAL  PERIOD  AT  RE 


START  PROGRAM 
PROGRAM  ORBMAINT 


( 


READ* , ALPHA 

ALPHAR=ALPHA*  PI / ISO . 0D+0 

I 

PRINT*  ,  '  ENTER  BAND,  •.  KM )  ‘ 

READ* , BAND 
3AND=BAND* 1000 . CD-t-G 


PRINT*, 'ENTER  THFAC ' 

READ* , THFAC 
THFAC=THFAC* 1 . 0D+0 

INITIALIZE  COUNTERS  USED  IN  PROGRAM 

;:,dex=o 

KOUNT= i 
TRAD=0 . 0D+0 


INITIALIZE  NONDIMENSIONALIZED  STATE  VARIABLES 


C 

C 

c 


c 

c 

c 


c 


c 

c 

c 


c 


c 

c 


XBAR ( 1 ) =R0/RE 
XBAR ( 2 ) =0D+0 
XBAR ( 3 ) =GD+0 

XBAR ( 4 ) =2 . GD+0  *PI *TAU/TAUR0 
COMPUTE  INITIAL  ORBITS  VELOCITY  AND  ENERGY 


V0= ( MU / RO )  *  *0 . 5 
E0=( (VO* VO ) / 2 . -MU/RQ ) 

COMPUTE  RADIUS, VELOCITY, ENERGY  FOR  DESIRED  BAND  MIN  AND  MAX 


RMIN=R0-BAND/2 . 0 

VMIN= ( MU/RMIN ) **0 . 5 

EMIN=( (VMIN*VMIN)/2. -MU/RMIN) 

RMAX=R0+BAND/2.0 

VMAX= ( MU / RMAX ) **0.5 

EMAX= ( (VMAX*VMAX)/2.0-MU/RMAX) 


SET  UP  OUTPUT  FILE  HEADERS 


I8AND=INT ( BAND ) 

I THMAX= I NT ( THMAX ) 
I ALPHA= I NT ( ALPHA ) 


WRITE ( 20 , * ) ' ! ! BAND= ' , IBAND , ' 
WRITE ( 20, *) ' ! ! EMIN= ' , EMIN, ' 
WRITE ( 30 , * ) ' ! ! BAND= ' , IBAND , ' 
WRITE ( 40, * ) ' ! !BAND=' , IBAND, ' 


THMAX= ' , ITHMAX, '  ALPHA= ' , I ALPHA 
EMAX= ' , EMAX 

THMAX=' , ITHMAX, '  ALPHA= ' , IALPHA 
THMAX= ' , ITHMAX , '  ALPHA= ' , IALPHA 


WRITE ( 20, * )  '  !  ! 

*  ANGM 

ORBITS  RADIUS 

ENERGY ' 

VELOCITY 

MFB 

MFBK 

WRITE (20, *  )  '  !  ! 

(KM) 

(KM/SEC) 

(KG) 

(KG)  ' 

WRITE ( 30, *)  '  !  ! 

*  PERIOD ' 

ORBITS  SMA 

ECC 

APOGEE 

PERIGEE 

WRITE ( 40, *)  ’  !  ! 
*  GAMMA' 

TIME  ORBITS 

DRAG 

THRUST 

MASS 
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o n  o  oooooo  noo  oonn  non 


C  CALCULATIONS 

C  CALL  DRAG  TO  SET  INITIAL  VALUE  FOR  FORCED  KEPLER I AN  MOTION 

C  USED  TO  COMPUTE  KEPLERIAN  THRUST 

CALL  Z RAG ( B  B AR , RE , X BAR , BETA , MB AR , V BAR , D ) 

DK=D 

TK=DK*MBAR*ME*RE/ ( TAU*TAU ) 

C 

c  CALL  SUBROUTINE  DRAG  TO  COMPUTE  PRESENT  DRAG 

100  CALL  DRAG  (  3BAR ,  RE , X8AR , BETA , MBAR , VBAR , D ) 

C 

C  CALCULATE  E  TO  SEND  TO  THRUST  SUBROUTINE  TO  DETERMINE 

C  IF  THRUSTER  SHOULD  BE  ON  OR  OFF 

C 

E= ( RE*VBAR/TAU ) **2 . 0/2 . 0D+0-MU/ (XBAR ( 1 ) *RE ) 

C 

C  CALL  SUBROUTINE  THRUST  TO  DETERMINE  IF  THRUSTERS  ARE  ON  OR 

C  OFF  AND  THE  SET  THE  VALUE  OF  THE  THRUST  TERM 

C 

CALL  THRUST ( E0 , XBAR , TH BAR , THFAC , RMIN , EMAX , E , ETH , THBM ) 

CALL  SUBROUTINE  EQN  TO  UPDATE  ORBITAL  EQUATIONS 

CALL  EQN ( D , XBAR , XBDOT , VBAR , THBAR , ALPHAR , MBAR , P I ) 

CALL  SUBROUTINE  RK4  TO  USE  RUNGE-KUTTA  FOR  INTEGRATION 
OVER  TIME 

CALL  RK4 ( TB , XBAR , XBDOT , J , TINCB , INDEX ) 

CHECK  TO  SEE  IF  RUNGE-KUTTA  COMPLETE 
IF (INDEX  .NE.  0)  GO  TO  100 


RUNGE-KUTTA  DONE  FOR  THIS  ITTERATION  OF  ORBIT, 

UPDATE  VARIABLES  FOR  CALCULATION  OF  ORBITAL 
PARAMETERS  AND  FOR  OUTPUT 

TH=THBAR*M£*RE/ ( TAU  *  TAU ) 

R=XBAR(1)*RE 
V=RE  * VBAR / TAU 

MF=TH*TINC/(SPI*G)  • 

M=M-MF 
MFT=MFT+MF 

GAMMAR=ATAN ( XBAR (2) / ( XBAR ( 1 ) *XBAR ( 4 ) ) ) 

GAMMAD*GAMMAR*180 . 0D+0/PI 
ANGM=R*V*COS ( GAMMAR ) 

ENERGY=< (V*V)/2.Q) - (MU/R ) 

TRAD»TRAD+TINCB  ft 

OANG=TRAD*2 . *PI 

COMPUTE  KEPLERIAN  FUEL  VALUES 

MFK=TK*TINC/ ( SPI *G ) 

MFTK-MFTK+MFK 

C  ft 


42 


CALL  SUBROUTINE  ORB PAR  TO  COMPUTE  ORBITAL  PARAMETERS 


IALL  ORBPAR < ENERGY , ANGM , R , MU ,  PI , ECO . SMA , APOGEE , 


COMPUTE  NEW  M3AR  AND  PUT  ?. ,  7  IN  OUTPUT  FORM 


MBAR=M/ME 
?.=R/ 1000.0 
7=V/ 1000.0 


DATA  OUT 


C  CHECK  TO  SEE  IF  TIME  TO  PRINT  OUTPUT 

IF ( KOUNT  .LT.  DNINT { PTIB/TINCB ) )  GO  TO  2C0 
C  PRINT  OUTPUTS  TO  OUTPUT  FILES  OUT,  RAT,  ORB 

c 

WRI TE ( 2  0 , 2 ) TB , R , V , MFT , MFTK , ANGM , E 
1  FORMAT  ( 2X , F6 . 2 , 2X , F9 . 3 , 2X , F7 . 4 , 2X, F8 . 2 , 2X , F8 . 2 , 2X , F12 . C , 

*2X , F10 . 0 ) 

WRITE  (  3 0 , 3  )  TB ,  SxMA ,  ECC ,  APOGEE ,  PERIGE ,  PERIOD 

3  FORMAT (2X,F7.2,iX,F10.3,3X,F4.3,3X,F10.3,iX,F10.3,2X,F8.4 > 


C 

4 

C 

C 

C 

C 

c 

c 

c 

c 

200 

c 

r 

C 

C 

C 

C 


WRITE  (  40 , 4  )  T ,  TB  ,  D  ,  TH ,  M ,  GAMMAD 

FORMAT (2X, F8.0, IX, F6.2 , 2X , F12 . 9 , 3X , F5 . 0 , 5X , F9 . 3 , IX, F10. 4) 
RESET  COUNTER 
KOUNT^O 

CHECK  TO  SEE  IF  ORBIT  COMPLETE,  RESET  ORBITAL  ANGLE 

IF  (TRAD  .GE.  1.0D+0)  TRAD=0.0D+0 

UPDATE  COUNTER 
KOUNT=KOUNT+l 

CHECK  FOR  END  OF  PROGRAM 

IF (TB  .LT.  TFB )  GO  TO  100 

END 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C 

c 

c 

c 


c 

c 

c 

c 

c 


THIS  SUBROUTINE  CALCULATES  DRAG  FOR  USE  IN  THE  EQUATIONS 
OF  MOTION. 

SUBROUTINE  DRAG ( BBAR, RE , XBAR, BETA, MBAR, VBAR, D ) 

IMPLICIT  REAL *8 ( A-H, M-Z ) 

DIMENSION  XBAR ( 4 ) 

SET  REFERENCE  ALTITUDE  FOR  ATMOSPHERIC  DENSITY  CALCULATIONS 
RR£F*=6.638145D+6 


43 


o  o  o  n  o  o  o  o  on  no  no  o  o  o  o  o  o  o  o  no 


COMPUTE  PRESENT  NONDIMENSIONALIZED  VELOCITY 
VBAR=<  ( XBAR  (  2  )  *  XBAR  (  2  )  )  +  ( XBAR<  1 )  *XBAR  (  4  )  }**2)”0.E 

COMPUTE  EXPONENTIAL  TERM  IN  DRAG  EQN 

Z=EXP I -BETA* RE* ( XBAR ( 1 ) - <  RREF/RE ) ) ) 

COMPUTE  DRAG  TERM 

D=E* VBAR* VBAR / ( 2 . 0*MBAR*BBAR ) 

RETURN 

END 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

THIS  SUBROUTINE  UPDATES  THE  EQUATIONS 
OF  MOTION 

SUBROUTINE  EQN ( D , XBAR , XBDOT , VBAR , THBAR , ALPHAR , MBAR ,?Z) 
IMPLICIT  REAL *8 ( A-H,M-Z ) 

DIMENSION  X3AR f  4 )  , XBDOT ( 4 ) 


A=XBAR ( 1 ) *XBAR ( 4 ) *XBAR ( 4 ) 
B=4.0*PI*PI/ (X8AR< 1) *XBAR(1) ) 
C=D  * { XBAR ( 2 ) / VBAR ) 
E=THBAR*SIN( ALPHAR) /MBAR 


F=2.0*XBAR(4)*XBAR(2)/XBAR(1) 
G=(D/XBARP ) ) *(XBAR(1)*XBAR(4)/VBAR) 
H=THBAR*COS (ALPHAR) / (MBAR*X3AR< 1 ) ) 


XBDOT ( 1 ) =XBAR ( 2 ) 
XBDOT ( 2 ) =A-B-C+E 
XBDOT ( 3 ) =X8AR ( 4 ) 
XBDOT ( 4 ) =-F-G+H 


RETURN 

END 


:CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

THIS  SUBROUTINE  CALCULATES  ORBITAL  PARAMETERS 

SUBROUTINE  ORBPAR( ENERGY , ANGM, R, MU , PI , ECC , SMA , APOGEE , PERIGE , 
* PERIOD) 

IMPLICIT  REAL  *  8 ( A-H , M- Z ) 

C  PROBLEM  COULD  ARISE  IN  COMPUTING  ECCENTRICITY  IF  VERY 

C  SMALL-  COULD  TRY  TO  TAKE  SQF.T  OF  NEG  NUMBER.  TO 

C  PREVENT,  IF  ECCENTRICITY  IS  LESS  THAN  IE-6 ,  SET 

C  EQUAL  TO  0. 

C 


CJUtN  rv)  oun  ro  UO 


?ROB  =  < 1. r2 . 'ENERGY *ANGM*ANGM/ (MU*MU)  ) 
IF  ( DABS ( PROB )  .  LT.  l.D-12)  THEN 

ECC=0 . CD+O 
ELSE 

ECC=?ROB*  *0 . 5 

ENDIF 


5MA= i -MU/ (2 . 'ENERGY \ ) / 1000.0 


SMAM=SMA* 1000.0 
APOG£E=SMA*( 1.0+ECC) 

PERIGE=SMA* ( 1 . 0-ECC 1 

PERIOD* 2 . 0*PI » ( SMAM"  1 . 5 ) / ( MU* *0 . 5 ) / 60 . 0 

Q 

RETURN 

END 

C 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C  THIS  SUBROUTINE  DOES  A  FOURTH  LEVEL  RUNGE-KUTTA 

-  INTEGRATION. 

C 

SUBROUTINE  RK4 < TB , XBAR , XBDOT , J , TINCB , INDEX) 

IMPLICIT  PEAL *8 ( A-H,M-Z ) 

INTEGER  INDEX, I 
C 

DIMENSION  XBAR ( 4 ) , XBDOT ( 4 ) , SAVED ( 4 ) , SAVEX ( 4 ) 

C 

INDEX=INDEX+I 
GO  TO  (1,2, 3, 4), INDEX 
C 

1  DO  10  1*1, J 

SAVEX  ( I )  =XBAP.  ( I ) 

SAVED ( I ) =XBDOT ( I ) 

10  XBAR ( I ) = SAVEX ( I ) +0 . SD+0*TINCB*XBDOT ( I ) 

TB=TB+0 . 5D+0*TINCB 
RETURN 


» 


» 


» 


DO  20  1*1, J 

SAVED ( I ) = SAVED ( I ) +2 . D+0*XBDOT ( I ) 

0  XBAR ( I ) =SAVEX( I ) +0. 5D+0*TINCB*XBDOT ( I ) 

RETURN 


DO  30  1=1, J 

SAVED ( I ) = SAVED ( I ) +2 . D+0*X8DOT ( I ) 
3  XBAR ( I ) =SAVEX ( I ) +TINCB*XBDOT ( I ) 

TB=TB+0. 5D+0*TINCB 
RETURN 


DO  40  1=1, J 

3  XBAR( I ) =SAVEX ( I ) +TINCB/ 6 . D+0* ( SAVED ( I ) +XBDOT ( I ) ) 

INDEX=0 
RETURN 


» 


» 
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on  nnpnnn  ooo  non  non  non  non 


END 


C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

THIS  SUBROUTINE  DETERMINES  IF  THRUSTERS  SHOULD  3E  ON  OR 
C  OFF  TO  TRY  TO  STAY  WITHIN  THE  DESIRED  BAND 

SUBROUTINE  THRUST ( EO , XBAR , THBAR , THFAC , RMIN, EMAX , E , ETH , 

■'THBM) 

IMPLICIT  REAL *8 ( A-H, M-Z ) 

r* 

DIMENSION  XBAR ( 4 }  , XBDOT ( 4 ) 

p 

RE=6. 3782D+6 

C  COMPUTE  NOND IMENS IONALI ZED  POINT  FOR  THRUSTERS  TO  TURN  ON 

C  BASED  ON  RADIUS. 

c 

RTH=RMIN/RE 

C 

C  ETH=EMAX- ( EMAX-EMIN ) *THFAC 

c 

IF  THRUSTERS  ARE  ON  GO  TO  SEE  IF  TIME  TO  TURN  OFF 
IF  (THBAR  .EQ.  THBM)  GO  TO  100 

CHECK  IF  RADIUS  BELOW  BAND 
IF  ( XBAR ( 1 )  .LE.  RTH )  THEN 

CHECK  IF  ENERGY  BELOW  INITIAL  ENERGY 
IF  (E  .LE.  EO)  THEN 

TURN  THRUSTERS  ON 
THBAR=THBM 

ELSE 

OTHERWISE  LEAVE  OFF 
THBAR=0D+0 

ENDIF 

RADIUS  NOT  BELOW  BAND,  LEAVE  OFF  "  RUSTERS 

ELSE 

THBAR=QD+0 

ENDIF 

IF  THRUSTERS  ARE  ON,  SEE  IF  TIME  TO  TURN  OFF 
)0  IF  (E  .GE.  EO)  THBAR=0D+0 


RETURN 

END 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


» 


» 
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» 


3. 


SAMPLE  INPUT  FILE 


o . 63ol4b0D+06 
2 . 00000000+04 
3 . OOOOOOOD+02 
0 . OQOOOOOD+OO 
1. OOOOOOOD+02 
2.0000000D-04 
1. OOOOOOOD-Ol 
2 . 5000000D+04 
2 . 12  OOOOOD-05 
3 . Q000000D+02 


C.  SAMPLE  OUTPUT  FILE 


:band  = 

25000 

THMAX= 

300 

!EMIN  = 

-30080181.48670808 

EMAX  = 

:  ORBITS 

RADIUS 

VELOCITY 

MFB 

(KM) 

(KM/SEC) 

(KG) 

1.00 

6636.385 

7.7501 

0.00 

2.00 

6634.580 

7.7512 

0.00 

3.00 

6632.737 

7.7523 

0.00 

4.00 

6630.855 

7.7534 

0.00 

5.00 

6628.930 

7.7545 

0.00 

6.00 

6626.951 

7.7557 

0.00 

7.00 

6634.414 

7.7513 

110.11 

8.00 

6628.660 

7.7600 

187.56 

9.00 

6625.232 

7.7641 

224.78 

10.00 

6622.517 

7.7673 

259.52 

11.00 

6620.638 

7.7695 

293.54 

12.00 

6619.667 

7.7706 

327.04 

13.00 

6619.631 

7.7706 

359.91 

14.00 

6620.517 

7.7696 

392.48 

15.00 

6622.261 

7.7676 

425.88 

16.00 

6624.814 

7.7642 

452.25 

17.00 

6628.234 

7.7599 

483.47 

18.00 

6632.050 

7.7555 

519.87 

19.00 

6636.064 

7.7509 

556.16 

20.00 

6640.091 

7.7462 

592.24 

21.00 

6643.950 

7.7417 

628.33 

22.00 

6647.468 

7.7377 

663.58 

23.00 

6650.501 

7.7342 

699.25 

24.00 

6652.923 

7.7314 

734.93 

ALPHA  =  "0 

-29967109.06483507 

MFBK  ANGM  ENERGY 

(KG) 

6.49  51432570538.  -30031067. 
12.98  51425887643.  -30038873. 
19.47  51418951743.  -30046977. 
25.96  51411742642.  -30055404. 
32.45  51404237633.  -30064181. 
38.94  51396411113.  -30073338. 
45.43  51425131788.  -30039612. 
51.91  51438552722.  -30023833. 
58.40  51438811678.  -30023532. 
64.89  51438814899.  -30023533. 
71.38  51438812800.  -30023539. 
77.87  51438816720.  -30023537. 
84.36  51438831447.  -30023522. 
90.85  51438826115.  -30023529. 
97.34  51438818018.  -30023539. 
103.83  51435956478.  -30026876. 
110.32  51434560951.  -30028501. 
116.81  51434837100.  -30028179. 
123.30  51435080714.  -30027896. 
129.79  51435311164.  -30027628. 
136.28  51435512393.  -30027394. 
142.76  51435699197.  -30027178. 
149.25  51435880508.  -30026968. 
155.74  51436047241.  -30026776. 


35.00 

no54.t>39 

"."'294 

~"0.60 

162.23 

51436208297 

;002ts59u. 

36.00 

0655.596 

".7283 

806.2" 

168.72 

51436378754. 

’0026392 

33700 

0655. r51 

".7281 

84 1 .63 

175.2: 

5 1436530808. 

38.00 

o655. 108 

".7289 

876.99 

181.70 

51436685225. 

'U026U38 

39.00 

o653.71 1 

".7306 

Jil.S3 

188.19 

5143685095b. 

.'002584b 

30.00 

6651.615 

"’.7330 

°47. 19 

194.68 

51437006143. 

-  .'0025656 

31.00 

6648.943 

".7362 

981.93 

201.17 

51437196213. 

;002544b 

32.00 

6645.784 

".7399 

1016.57 

207.66 

51437368920. 

-30025245 

33.00 

6642.315 

".7440 

1051.21 

214.15 

51437587533. 

-.'0024991 

34.00 

6638.654 

".7483 

1085.74 

220.64 

51437799569. 

-30024744 

35.00 

6634  978 

7.7526 

1120.38 

227.13 

51438037086. 

-30024457 

36.00 

6631.441 

".7568 

1154.81 

233.62 

51438289123. 

-30024173 

37.00 

6628.199 

7.7606 

1189.65 

24C.10 

51438574149. 

-30023 840 

38.00 

6625.379 

7.7639 

1223.36 

246.59 

51438836529. 

-30023533 

39.00 

6623.118 

"77666 

1255.62 

253.08 

51438840830. 

30023530 

40.00 

6621.515 

".7684 

1286.84 

259.57 

51438834154. 

-30023540. 

41.00 

6620.631 

7.7695 

1317.86 

266.06 

51438836523. 

-30023538. 

42.00 

6620.505 

7.7696 

1348.26 

272.55 

51438836153. 

-3°023539. 

43.00 

6621.127 

7.7689 

1379.38 

279.04 

51438840521. 

-30023533. 

44.00 

6622.477 

7.7673 

1410.30 

285.53 

51438837255. 

-30023536. 

45.00 

6624.516 

7.7646 

1436.45 

292.02 

51436551488. 

-30026202. 

46.00 

6627.376 

7.7609 

1462.20 

298.51 

51434160812. 

-30028987. 

47.00 

6630.657 

7.7571 

1496.42 

305.00 

51434443176. 

-30028656. 

48.00  6634.186  7.7530  1529.72  311.49  51434694318.  -30028363. 

49.00  6637.818  7.7488  1563.32  317.98  51434926906.  -30028091. 

50.00  6641.411  7.7446  1597.65  324.47  51435148627.  -30027833. 

51.00  6644.824  7.7407  1631.46  330.95  51435348781.  -30027600. 

52.00  6647.925  7.7371  1664.65  337.44  51435535867.  -30027382. 

53.00  6650.596  7.7340  1698.46  343.93  51435712082.  -30027177. 

54.00  6652.744  7.7315  1732.27  350.42  51435892956.  -30026967. 

55.00  6654.278  7.7298  1765.66  356.91  51436042361.  -30026794. 

56.00  6655.151  7.7288  1799.68  363.40  51436201010.  -30026610. 

57.00  6655.343  7.7286  1833.38  369.89  51436371724.  -30026411. 

58.00  6654.841  7.7292  1867.40  376.38  51436530981.  -30026227. 

59.00  6653.669  7.7306  1901.31  382.87  51436693930.  -30026037. 

60.00  6651.866  7.7327  1934.71  389.36  51436843665.  -30025863. 

61.00  6649.535  7.7354  1969.04  395.85  51437033997.  -30025642. 

62.00  6646.736  7.7387  2003.05  402.34  51437203180.  -30025445. 

63.00  6643.614  7.7424  2036.66  408.83  51437417724.  -30025195. 

64.00  6640.256  7.7464  2070.47  415.32  51437615550.  -30024965. 

65.00  6636.822  7.7504  2103.76  421.80  51437843332.  -30024699. 

66.00  6633.438  7.7544  2136.74  428.29  51438096137.  -30024404. 

67.00  6630.226  7.7582  2170.24  434.78  51438361782.  -30024093. 


48 


b8.uO 

ot>27. 315 

"."615 

2203.95 

441.27 

51438659553. 

;002 3 "4 5. 

-v9.00 

'-624. 808 

"."645 

2236.21 

-147. "6 

5 1438839896. 

•  >025554 

"O.UO 

^o22.852 

"."669 

2267.43 

454.25 

'1438841 ,4" 

-•02553- 

"1.00 

••o2 1 .455 

"."685 

2298.  ;4 

-60.74 

5U3885405V 

'9025520. 

"z.oo 

i'b20."30 

"."694 

2328.33 

-67.23 

5 1438848302 

•9023528. 

"'.00 

n620.676 

"."694 

2358.32 

473. "2 

5 1438838835. 

.'0023539. 

"4.00 

6621.291 

".7687 

2389.23 

480.21 

51438840540 

902353" 

"5.00 

0622.55 1 

".7672 

2419.84 

486.70 

51438841940. 

-9023534 

"6.00 

0624.435 

".7647 

2446.  iO 

493.19 

51436656305. 

390260" ; . 

"7.00 

6627.090 

".7612 

2471.84 

499.68 

51433850507. 

-30029352. 

"8.00 

6630.147 

".7576 

2505.45 

506.17 

51434145633. 

-30029006. 

"9.00 

5633.448 

".7538 

2538.84 

512.66 

51434396247 

.'002 S":  3. 

80.00 

6636.870 

"  7499 

2572.24 

519.14 

5 14346425 "4. 

-30028425. 

81.00 

6640.288 

7.7459 

2605.53 

525.63 

51434851838. 

■30028:81. 

32.00 

6643.582 

".7421 

2639.14 

532.12 

51435066093. 

-  '002793 i . 

83.00 

5646.639 

7.7386 

2672.94 

538.51 

51 4352604 i 9. 

3002~"05. 

84.00 

5649.350 

".7354 

2706.03 

545.10 

51435436862. 

-50027500. 

85.00 

6651.626 

".7328 

2739.84 

551.59 

51435610960. 

-'002"29“ 

86.00 

6653.393 

".7308 

2773.55 

558.08 

51435778411. 

-50027103. 

87.00 

6654.598 

".7294 

2806.74 

564.57 

51435948364. 

-50026905. 

88.00 

6655.195 

7.7287 

2840.13 

571.06 

51436102292. 

-3002672". 

89.00 

6655. 170 

7.7288 

2874.05 

577.55 

51436255865. 

■30026549. 

90.00 

6654.533 

7.7295 

2907.24 

584.04 

51436418385. 

-30026360. 

91.00 

6653.305 

7.7310 

2940.53 

590.53 

51436578220. 

-3002  61 74. 

92.00 

6651.546 

7.7331 

2974.03 

597.02 

51436757933. 

-50025965. 

93.00 

6649.293 

7.7357 

3007.74 

603.51 

51436918736. 

-30025778. 

94.00 

6646.656 

7.7388 

3041.55 

609.99 

51437110521. 

-30025555. 

95.00 

6643.718 

7.7423 

3075.35 

616.48 

51437314954. 

-30025317. 

96.00 

6640.575 

7.7460 

3109.16 

622.97 

51437527062. 

-30025070. 

97.00 

6637.340 

7.7498 

3142.25 

629.46 

51437748538. 

-30024811. 

98.00 

6634.132 

7.7536 

3175.44 

635.95 

51437987351. 

-50024532. 

99.00 

6631.065 

7.7572 

3208.84 

642.44 

51438257984. 

-30024216. 

100.00 

6628.231 

7.7605 

3242.54 

648.93 

51438533304. 

-30023894 
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APPENDIX  B 


Figure  4 
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Figure  8 
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Figure  4 


» 
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APPEKDIX  D 


Figure  1  a=60° 
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Figure  2  a=65° 
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Figure  4  a=60° 
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APPENDIX  F 


Figure  2 
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Figure  3 
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APPENDIX  G 


Figure  1 


Figure  2 
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Figure  3 
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APPENDIX  H 


Figure  1 


Figure  2 
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